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Abstract  An  advanced  model  of  irreversible  thermoelectric  generator  with  a  generalized  heat  transfer 
law  is  established  based  on  finite  time  thermodynamics.  The  generalized  heat  transfer  law  represents  a 
class  of  heat  transfer  laws  including  Newtonian  heat  transfer  law,  linear  phenomenological  heat  transfer 
law,  radiative  heat  transfer  law,  Dulong-Petit  heat  transfer  law,  generalized  convective  heat  transfer  law 
and  generalized  radiative  heat  transfer  law.  The  inner  effects  including  Seebeck  effect,  Fourier  effect, 
Joule  effect  and  Thomson  effect,  and  external  heat  transfer  are  taken  into  account  in  the  model.  The 
Euler-Lagrange  functions  at  maximum  power  output  and  maximum  efficiency  are  established.  Applying 
the  model  to  a  practical  example  in  engineering,  it  is  found  that  the  external  heat  transfer  law  does  affect 
the  characteristics  and  optimal  performance  of  the  thermoelectric  device,  and  the  maximum  power  output 
and  maximum  efficiency  with  Newtonian  heat  transfer  law  are  the  maximum  among  the  several  typical 
heat  transfer  laws.  The  results  can  offer  principles  for  the  power  and  efficiency  optimization  of  practical 
thermoelectric  generators  at  various  external  heat  transfer  conditions. 
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1.  Introduction 

The  recent  worldwide  energy  consciousness  has  focused 
attention  on  unconventional  means  of  producing  electrical 
energy.  The  two  of  these  -  the  photovoltaic  process  and  the 
thermoelectric  process  -  are  the  only  processes  that  convert 
solar  energy  directly  into  electrical  energy  [1-4].  Thermoelec¬ 
tric  generation  technology,  due  to  its  several  kinds  of  merits, 
especially  its  promising  applications  to  waste  heat  recovery,  is 
becoming  a  noticeable  research  direction  [5-7].  Considerable 
efforts  have  been  developed  to  identify  new  classes  of  ther¬ 
moelectric  materials.  In  addition  to  the  improvement  of  ther¬ 
moelectric  materials,  the  system  modeling  and  optimization  of 
thermoelectric  generators  are  equally  important  in  improving 
the  performance  of  thermoelectric  generators  [8-10]. 
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So  far,  most  of  the  researches  on  thermoelectric  gener¬ 
ators  are  mainly  based  on  the  non-equilibrium  thermody¬ 
namics  [11-14]  and  the  finite  time  thermodynamics  (FTT) 
[15-22].  Researches  based  on  the  non-equilibrium  thermody¬ 
namics  paid  attention  to  the  inner  structure  of  thermoelectric 
elements.  A  significant  increase  in  the  power  output  from  a 
module  can  be  achieved  by  modifying  the  size  of  the  thermo¬ 
electric  elements  [23].  Burshtein  [24],  Chen  et  al.  [25],  Omer 
and  Infield  [26],  Moizhes  et  al.  [27,28],  Seifert  et  al.  [29,30], 
Drabkin  and  Ershova  [31]  and  other  researchers  investigated 
the  power  output  or  efficiency  of  single-element  thermoelectric 
generator.  Actually,  a  commercial  thermoelectric  generator  is  a 
multi-element  device,  which  is  composed  of  many  fundamental 
thermoelectric  elements.  Rowe  [32],  Sisman  and  Yavuz  [33], 
Mayergoyz  and  Andrei  [34]  investigated  the  power  output  and 
efficiency  of  a  multi-element  thermoelectric  generator  in  vari¬ 
ous  working  conditions.  However,  the  external  heat  transfer  ir¬ 
reversibility  was  not  considered  in  the  literature  hereinbefore. 

Thermoelectric  devices  should  be  connected  with  heat 
exchangers  to  absorb  and  dissipate  heat  [1-10].  Much  work  has 
shown  that  the  heat  transfer  irreversibility  between  an  engine 
and  its  external  reservoirs  affects  the  thermodynamic  processes 
obviously.  The  theory  of  finite  time  thermodynamics  is  a 
powerful  tool  for  the  performance  analysis  of  thermodynamic 
processes.  Some  authors  have  investigated  the  performances 
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Nomenclature 

A  Area  (m2) 

C  Coefficient 

ej  Electrical  current  density  (A/m2 ) 

/  Ratio  of  thermal  conductance  allocation 

I  Working  electrical  current  (A) 

I<  Thermal  conductance  ( WK“ 1 ) 

k  Thermal  conductivity  ( Wm_  1 1<_  1 ) 

L  Length  (m) 

N  Number 

P  Power  output  (W) 

Q.  Heat  flow  rate  (W) 

R  Electrical  resistance  (Q) 

T  Temperature  (K) 

Greek  letters 

a  Seebeck  coefficient  (VK-1 ) 

A  Difference 

/x  Thomson  coefficient  (VK-2) 

p  Electrical  resistivity  (£2m) 

a  Electrical  conductivity  (£2_1m_1) 

Subscripts 

c  Cold  junction 

H  Heat  source 

h  Hot  junction 

L  Heat  sink 

n  N-type  semiconductor  leg 

p  P-type  semiconductor  leg 

P  At  maximum  power 

T  Total 

rj  At  maximum  efficiency 


of  thermoelectric  generators  using  a  combination  of  finite 
time  thermodynamics  and  non-equilibrium  thermodynamics. 
Gordon  [35,36],  Yan  and  Chen  [37],  Agrawal  and  Menon  [38], 
and  Chen  and  Wu  [39]  analyzed  the  performance  of  single¬ 
element  thermoelectric  generators  with  finite  rate  heat  transfer 
between  the  thermoelectric  device  and  its  external  heat 
reservoirs.  Further,  Chen  et  al.  [40],  Yu  and  Zhao  [41],  Hsiao 
et  al.  [42],  Astrain  et  al.  [43],  and  Meng  et  al.  [44]  investigated 
the  characteristics  of  multi-element  thermoelectric  generation 
system  with  external  heat  transfer. 

In  the  analyses  mentioned  above,  the  heat  transfer  between 
the  thermoelectric  generator  and  the  heat  reservoirs  was  as¬ 
sumed  to  obey  Newtonian  (linear)  heat  transfer  law  Q.  oc  AT. 
Much  work  has  shown  that  the  heat  transfer  law  between  the 
device  and  its  external  reservoirs  affects  the  performances  of 
thermodynamic  cycles  strongly.  However,  when  the  external 
heat  reservoirs  obey  nonlinear  heat  transfer  laws,  the  analytical 
results  are  hard  to  obtain.  Chen  et  al.  [45]  and  Meng  et  al.  [46] 
established  a  finite  time  thermodynamic  model  of  a  multi¬ 
element  thermoelectric  generator  with  linear  phenomenolog¬ 
ical  law  and  radiative  law  and,  consequently,  investigated  the 
characteristics  of  the  device. 

Reviewing  the  former  literature  concerning  thermoelectric 
generators,  some  features  can  be  concluded  as  follows: 

(1)  Researches  concerning  the  inner  structure  of  a  thermoelec¬ 
tric  generator  are  without  considering  the  external  heat 

transfer  irreversibility. 


(2)  Researches  concerning  the  external  heat  transfer  optimiza¬ 
tion  of  a  thermoelectric  generator  are  without  considering 
the  temperature  dependence  of  thermoelectric  properties, 
i.e.  the  Thomson  effect  of  thermoelectric  materials. 

(3)  Most  of  the  researches  assumed  that  the  external  heat 
transfer  obeys  Newtonian  heat  transfer  law  but  not 
concerned  nonlinear  heat  transfer  law. 

The  heat  transfer  law  represents  the  characteristic  and 
regularity  of  the  transfer.  The  laws  would  be  complex  and  the 
well-known  Newtonian  heat  transfer  law  is  not  the  only  law  in 
practical  heat  transfer  process.  There  are  many  other  nonlinear 
heat  transfer  laws  and  the  heat  transfer  law  does  affect  the 
performance  of  the  device.  Many  researchers  have  concerned 
the  characteristics  of  conventional  heat  engines,  refrigerators 
and  heat  pumps  with  linear  phenomenological  heat  transfer 
law  Q.  oc  AT-1  [47,48],  radiative  heat  transfer  law  Q.  a  AT4 
[49,50],  Dulong-Petit  heat  transfer  law  Q.  oc  (AT)1-25  [51,52] 
generalized  convective  heat  transfer  law  Q  oc  (AT)m  [53,54] 
and  generalized  radiative  heat  transfer  law  Q.  oc  A  (Tn)  [55-57]. 
To  establish  a  general,  natural  expression  of  various  heat 
transfer  laws,  a  generalized  heat  transfer  law  Q.  oc  (A Tn)m, 
which  represents  a  class  of  heat  transfer  laws  including  the 
heat  transfer  laws  mentioned  above,  was  put  forward  by  Chen 
et  al.  [58]  and  extended  by  Li  et  al.  [59-61].  The  effects  of 
heat  transfer  laws  on  the  performance  characteristics  of  the 
thermodynamic  systems  were  investigated. 

However,  there  is  a  lack  of  a  model  of  thermoelectric  gen¬ 
erator,  which  can  be  applied  to  the  analysis  and  optimization 
of  a  multi-element  thermoelectric  generator  with  generalized 
external  heat  transfer  law,  and  can  be  applied  to  analyze  the  ef¬ 
fects  of  heat  transfer  laws  on  the  maximum  power  output,  max¬ 
imum  efficiency  and  the  optimal  variables.  This  paper  aims  to 
establish  an  advanced  model  of  an  irreversible  multi-element 
thermoelectric  generator  with  generalized  external  heat  trans¬ 
fer  law  Q.  a  (A Tn)m  by  introducing  the  theory  and  method 
used  in  literature  mentioned  above.  The  Fourier  heat  leakage 
loss,  Joule  heat  loss,  temperature  dependence  of  thermoelectric 
properties,  especially,  the  Thomson  heat  loss  are  taken  into  ac¬ 
count  in  the  model.  The  maximum  power  output  and  maximum 
efficiency  are  explored.  A  practical  example  in  engineering  is 
proposed  to  analyze  the  effects  of  several  typical  heat  transfer 
laws  on  the  optimal  variables,  i.e.  working  electrical  current  and 
ratio  of  thermal  conductance  allocation  of  the  heat  exchangers. 

2.  Thermodynamic  model  and  solution 

2.1.  An  advanced  finite  time  thermodynamic  model 

Figure  1  shows  the  advanced  finite  time  thermodynamic 
model  of  an  irreversible  thermoelectric  generator  with  external 
generalized  heat  transfer  law  established  in  this  paper. 
The  thermoelectric  generator  consists  of  P-type  and  N-type 
semiconductor  legs.  The  number  of  thermoelectric  elements 
is  N.  The  junctions  of  the  thermoelectric  elements  are  fixed 
at  a  thermal  conducting  and  electrical  insulting  ceramic  plate. 
The  thermoelectric  elements  are  insulated,  both  electrically  and 
thermally,  from  its  surroundings.  The  internal  irreversibility 
is  caused  by  Joule  loss  and  Fourier  heat  conduction  loss 
through  the  semiconductor  between  the  hot  and  cold  junctions. 
The  Joule  loss  generates  an  internal  heat  I2R}  where  R  is 
the  total  internal  electrical  resistance  of  the  thermoelectric 
element  and  I  is  the  electrical  current  generating  from  the 
thermoelectric  elements.  The  conduction  heat  loss  is  I<(Th  — 
Tc),  where  I<  is  the  thermal  conductance  of  the  thermoelectric 
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Figure  1 :  Finite  time  thermodynamic  model  of  an  irreversible  thermoelectric 
generator  with  generalized  heat  transfer  law. 

element,  Th  is  the  hot  junction  temperature  and  Tc  is  the  cold 
junction  temperature.  The  heat  flow  rates  through  the  hot  and 
cold  junctions  of  the  thermoelectric  elements  are  Q/,  and  Qc, 
respectively. 

The  external  irreversibility  is  caused  by  the  finite  rate  heat 
transfer  between  the  thermoelectric  generator  and  its  heat 
reservoirs.  The  temperatures  of  the  heat  source  and  heat  sink 
are  TH  and  TL ,  respectively.  The  thermal  conductance  of  the 
heat  exchangers  are  I<H  and  KL,  respectively.  The  total  thermal 
conductance  I(T  =  I<h  +  KL  of  the  external  heat  transfer  of 
the  thermoelectric  device  is  finite  and  fixed.  The  heat  flow  rate 
absorbed  from  the  heat  source  to  the  thermoelectric  generator 
is  Qh.  The  heat  flow  rate  dissipated  from  the  thermoelectric 
generator  to  the  heat  sink  is  Q.. 

As  a  commercial  thermoelectric  generation  module  is 
thermal  insulation  packaged,  the  heat  leakage  through  the 
lateral  face  can  be  neglected.  At  a  steady  work  state,  the 
temperature  distribution  of  the  air  gap  is  the  same  with 
the  thermoelectric  elements,  so  the  heat  transfer  of  the 
whole  device  can  be  treated  as  one  dimensional  heat  transfer 
approximately.  The  increment  rate  of  inner  energy  of  an 
infinitesimal  is  zero  at  steady-state,  so  one  can  obtain  the 
energy  conservation  equation  of  a  semiconductor  leg  as 
follows: 

Ql<in  ~  QlCout  +  Q/z  +  0/  =  0,  (1) 


coefficient,  cross  section  area,  length,  temperature  and  elec¬ 
trical  current  density  of  the  P-type  and  N-type  semiconductor 
leg.  Reforming  Eqs.  (2)  and  (3)  one  can  obtain  the  heat  conduc¬ 
tion  differential  equations  of  P-  and  N-type  semiconductor  legs 
(Eqs.  (4)  and  (5))  with  boundary  conditions  (Eqs.  (6)  and  (7))  as 
follows  [12]: 


dll 


dTr, 


V^+Mp/^f  +  —  -0, 


crpAp 


.  A  dT2  dTn  I2 
Mn  2  Mr>/  + 

U  A  U  A  O  n/\n 


o, 


(4) 

(5) 


TP(  0)  =  T„(0)  =  Tc,  (6) 

Tp(Lp)  =  Tn(Ln )  =  Th,  (7) 

where  kp,  op,  np,  Ap,  Lp,  Tp  and  kn,  on,  \in,  An,  Ln,  Tn  are  the 
thermal  conductivity,  electrical  conductivity,  Thomson  coeffi¬ 
cient,  cross  section  area,  length  and  temperature  of  the  P-type 
and  N-type  semiconductor  leg.  I  =  eJpAp  =  eJnAn  is  the  electri¬ 
cal  current.  In  Eqs.  (4)  and  (5),  the  first  term  (second  derivative) 
is  caused  by  thermal  conduction  of  P-  and  N-type  semiconduc¬ 
tor  legs ;  the  second  term  ( first  derivative )  is  caused  by  Thomson 
effect;  the  third  term  is  caused  by  Joule  effect.  The  Seebeck  ef¬ 
fect  is  the  surface  effect,  so  it  is  not  embodied  in  the  equations 
but  it  causes  the  electrical  current  through  the  elements. 

Taking  effect  of  temperature  dependence  of  thermoelectric 
properties  into  account,  kp,op  and  /zp  are  function  of  Tp\ 
while  kn,  on  and  /xn  are  function  of  Tn.  However,  such  a 
differential  equation  cannot  be  solved  analytically  in  general. 
Replacing  k,  a  and  /z  with  the  averaged  coefficients  k,  o  and  JL 
approximately  [62]  gives  approximations  of  Eqs.  (4)  and  (5)  as 
follows: 


—  dll 


din 
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(9) 


where /<p  ~  kp|r=(rh+rc)/2,  op  ~  op |r=(Th+Tc)/2,  and jiv  ~ 
/Zp|T=(Th+Tc)/2  for  P-type  semiconductor  legs  and  kn  = 
kn\T=(Th+Tc)/2,  on  =  <7n|T=(Th+Tc)/2,  and /zn  =  fin\T=(Th+Tc)/2  for 

N-type  semiconductor  legs.  The  total  heat  flow  rates  through 
the  hot  and  cold  junctions  are: 


where  Qi<in,  Qicout,  Q/x  and  Q/  are  the  Fourier  heat  input  (heat 
input  to  the  hot  junction  of  the  element  by  heat  conduction), 
Fourier  heat  output  (heat  output  from  the  cold  junction  of 
the  element  by  heat  conduction),  generated  Thomson  heat 
(generated  heat  because  of  Thomson  effect)  and  generated  Joule 
heat  (generated  heat  when  electrical  current  flows  through  the 
conductor).  Eq.  (1)  can  be  expressed  as: 

d  dT 

—  (Tp  +  dlp)kpAp  —  kpAp  +  fipeJpApdTp 
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On  In 

where  kp,  Op,  fip,  Ap,  Lp,  Tp,  cjp  and  kn,  on,  /zn,  An,  Ln,  Tn,  cjn 

are  the  thermal  conductivity,  electrical  conductivity,  Thomson 
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—  dTn 
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dx 
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T-  dTn 

+  knAn  — — 
dx 


(10) 


(11) 


where  ap  and  an  are  the  Seebeck  coefficients  of  the  P-  and 
N-type  semiconductor  legs,  and  the  subscript  h  and  c  represent 
the  hot  and  cold  sides.  In  Eqs.  (10)  and  (11),  the  first  term  is 
caused  by  Seebeck  effect,  and  the  second  and  the  third  terms  are 
caused  by  thermal  conduction  of  P-  and  N-type  semiconductor 
legs.  The  Joule  and  Thomson  effects  are  volume  effect,  so 
they  are  not  embodied  in  the  equations  but  they  affect  the 
temperature  distribution  of  the  elements. 
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It  is  assumed  that  the  heat  transfer  between  the  junctions 
of  the  thermoelectric  elements  and  their  respective  reservoirs 
obeys  a  generalized  heat  transfer  law  Q.  oc  (ATn)m  (including 
Newtonian  heat  transfer  law  Q.  oc  AT,  linear  phenomenological 
heat  transfer  law  Q.  oc  AT-1,  radiative  heat  transfer  law  Q  oc 
AT4,  Dulong-Petit  heat  transfer  law  Q.  oc  (AT)125,  generalized 
convective  heat  transfer  law  Q  oc  (AT)171  and  generalized  ra¬ 
diative  heat  transfer  law  Q.  oc  ATn),  as  did  the  finite  time  ther¬ 
modynamics  for  various  thermodynamic  processes  and  devices. 
According  to  the  theory  of  finite  time  thermodynamics,  one  has: 


To  reduce  the  cost  of  manufacture,  the  P-  and  N-type 
semiconductor  legs  are  made  with  same  sizes,  i.e.  Ap  =  An  =  A 
andLp  =  Ln  =  L.  So  one  has  knorn/(kporp )  =  1  by  Eq.(  18).  Similar 
doped  alloys  are  adopted  to  make  P-  and  N-type  semiconductor 
legs.  That  is  op  =  on  =  cr,  kp  =  kn  =  k,  ap  =  —an,  and 
l_ip  =  —\xn.  This  choice  is  only  due  to  the  simplification  of  the 
calculations.  According  to  Taylor’s  formula,  when  |x|  <$C  1,  ex  ~ 
1  +x  holds  true.  Based  on  above  assumptions,  Eqs.  (10)  and  (11) 
can  be  approximated  and  simplified  by  a  polynomial  expression 
of  current  /as  follows  at  the  case  of  JZI/K  <£.  1 


QH=KH(TZ-TZ)m  (12) 

QL  =  KL(Tcn-TLn)m.  (13) 

Eqs.  (12)  and  (13),  include  various  heat  transfer  laws.  Specially: 

(1)  If  n  =  1,  the  law  is  generalized  convective  heat  transfer 
law  [53,54].  Further, 

(a)  If  m  =  1,  the  law  is  Newtonian  heat  transfer  law. 

(b)  If  m  =  1.25,  the  law  is  Dulong-Petit  heat  transfer  law 
[51,52]. 

(2)  If  m  =  1,  the  law  is  generalized  radiative  heat  transfer 
law  [55-57].  Further, 

(a)  If  n  =  —1,  the  law  is  linear  phenomenological  heat 
transfer  law  [47,48]. 

(b)  If  n  =  4,  the  law  is  radiative  heat  transfer  law  [49,50]. 
The  energy-balance  equations,  i.e.  the  control  equations, 
are: 
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(14) 


(15) 


2.2.  Maximum  power  output  and  maximum  efficiency 


The  temperature  distributions  of  the  P-  and  N-type  semicon¬ 
ductor  legs  can  be  solved  by  Eqs.  (8)  and  (9)  as  follows: 

Tr(x)  =  T,-Fr,+  Tj0: (16) 

T,(«)  =  T,  +  F,x  +  ^  ~  y*  (e~  -  1).  (17) 

where  cop  —  [ipI/(kpApf  Fp  =  I/(<7p/xpApf  con  =  /xn//(/cnAn) 
and  Fn  =  I/(crnfinAn).  Th  and  Tc  are  not  fixed  because  they  can 
change  if  the  current  changes. 

For  maximum  figure  of  merit  of  the  thermoelectric  element 
Z  =  a2 the  size  of  the  thermoelectric  element  and  the 
physical  property  of  the  material  should  satisfy  the  following 
equation: 


^pLn  _  Kgn 

^ntf}  kv0v 


(18) 


Q,,  =  N[ahITh  +  K(Th  -  Tc) 

-0.5I2R-0.5JZI(Th-Tc)],  (19) 

Qc  =  N[acITc  +  K(Th  -  Tc) 

+  0.5 12R  +  0.5JZI(Th  -  Tc)]  (20) 

where  ah  =  cq,p  —  and  ac  =  acp  —  acn  are  the  Seebeck 
coefficients  of  the  thermoelectric  elements  at  hot  and  cold 
sides.  ~jl  =  JLp  —  JLn  =  2]Tf,  I<  =  I<p  +  I<n  and  R  =  Rp  +  Rn 
are  the  total  Thomson  coefficient,  thermal  conductance  and 
electrical  resistance  of  a  thermoelectric  element.  alT,  KAT,  I2R 
and  JLl  AT  are  the  rates  of  Peltier  heat,  Fourier  heat,  Joule  heat 
and  Thomson  heat,  respectively. 

I<  and  R  are  given  by: 


I(  =  I<p  +  Kn  =  kpAp/Lp  +  knAn/Ln  =  2kA/L,  (21) 

R  =  Rp+Rn=  Lp/{a~pAp )  +  Ln/(crnAn)  =  2  L/(WA).  (22) 

The  power  output  is  given  by: 

P  =  kh(TS  -  Thn)m  -  ~  Tf)m.  (23) 

The  efficiency  p  =  P/Qh  is  given  by: 


*7  =  1- 


I<L(Tcn 


Tf)m 


Tn)m- 


(24) 


The  power  output  of  a  thermoelectric  generator  depends  on 
the  working  electrical  current  /.  For  a  given  thermoelectric  gen¬ 
erator,  there  exists  an  optimal  electrical  current  correspond¬ 
ing  to  the  maximum  power  output,  which  has  been  proved  by 
much  literature.  Thus,  the  optimization  of  electrical  current  is  a 
basic  problem  of  the  internal  optimization  of  a  thermoelectric 
generator. 

More  external  thermal  conductance  means  better  heat 
exchange  but  bulkier  structures  and  higher  cost.  When  the 
external  heat  transfer  is  fixed,  the  allocation  of  thermal 
conductance  of  heat  exchangers  would  affect  the  performance 
of  the  device.  Therefore,  it  is  an  important  problem  that  how 
to  allocate  fixed  total  external  thermal  conductance  among 
hot  and  cold  sides  of  the  device  for  maximum  power  output 
or  maximum  efficiency.  The  optimum  allocation  of  thermal 
conductance  for  conventional  power  and  refrigeration  plants 
was  first  advanced  by  Bejan  [16,63,64]  and  advanced  by  many 
researchers.  To  describe  the  allocation,  a  ratio  of  thermal 
conductance  allocation  is  introduced  as/  =  I<H /(I<H  +I<L ).  Then 
one  has  I<H  =  fI<T  and  I(L  =  (1  -  f)I<T. 

Much  work  has  shown  that  when  the  external  heat  transfer 
law  is  linear  i.e.  Q  oc  AT,  there  will  exist  optimal  electrical 
currents  and  optimal  ratios  of  thermal  conductance  allocation 
corresponding  to  the  maximum  power  output  and  maximum 
efficiency,  respectively.  This  research  will  prove  that  when  the 
external  heat  transfer  law  is  nonlinear  i.e.  Q  oc  (A Tn)m  (n  / 
1  or  m  /  1),  there  will  also  exist  optimal  electrical  currents  and 
optimal  ratios  of  thermal  conductance  allocation. 
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The  Euler-Lagrange  function  for  maximum  power  output 
can  be  established  by  combining  Eq.  (23)  with  Eqs.  (14)  and  (15) 
as  follows: 

Lp  =  I <Tf(TS  -  K)m  -  I(T(  1  — /)(T"  -  T/T 

+  Xp(Krf(T"  -  Thn)m  -  N(ah  ITh  +  K(Th  -  Tc) 

-  0.5 12R- 0.5  JH(Th  —  Tc))) 

+  (/>p(I<tCI  ~f)(T"  -  T")m  -  N(acITc 
+  K(Th  -  Tc)  +  0.5 12R  +  0.5  JZI(Th  -  Tc))),  (25) 


-Ml  —  f  )mnT"~lDc(T"  -  T")-1)  =  0, 

(36) 

f~2DcDh 1  +  X,KTDh  +  <pnKTDc  =  0, 

(37) 

—XnN  (ah  I/,  —  IR  —  0.5  ]I(Th  —  Tc)) 

+  (prjN  (ac  TC+IR  +  0.5  JZ  (Th  -  Tc))  =  0, 

(38) 

where: 

°h  =  (t„  -  r„")m , 

(39) 

Dc  =  (t"  -  r")m . 

(40) 

where  kP  and  0P  are  the  Lagrangian  multipliers,  which  are 
constants  to  be  determined.  The  optimal  ratio  of  thermal 
conductance  allocation,  fP,  optimal  electrical  current,  IP,  and  the 
corresponding  hot  junction  temperature,  Th,  and  cold  junction 
temperature,  Tc,  should  satisfy  the  Euler-Lagrange  equations: 


3LP  3LP  dLp  3LP 

-  =  0,  - =  0,  -  =  0,  -  = 

3Th  3TC  3/  3/ 

Eq.  (26)  can  be  calculated  as: 

(1  +  kP)Ch  +  XPN  (ah  I  +  K  -  0.5 JZI) 

—  cj)PN  (K  +  0.5  JZI)  =  0, 

(l+0p)Cc  +  APN(-/C  +  O.5/I/) 

—  0PN  (ctc  /  -K-  0.5  JZI)  =  0, 

(1  +  XP)  (1*  -  TH)m  +  (1  +  <pP)  (T"  -  T")m  =  0, 
Xp  [Th  ah-IR-  0.5  JZ  (Th  -  Tc)] 

—  (pP  [Tc  <xc  +  IR  +  0.5  p  (Th  —  Tc)]  =  0, 

where: 


(26) 

(27) 

(28) 

(29) 

(30) 


Q,  =  KTfmnT”-'  (T"  -  T£)m  1 ,  (31) 

Q  =  Kt  (1  -/)  mnT"~J  (Tcn  -  .  (32) 

The  Euler-Lagrange  function  for  maximum  efficiency  can  be 
established  by  combining  Eq.  (24)  with  Eqs.  (14)  and  (15)  as 
follows: 


It  is  important  to  note  that  not  only  the  external  heat  trans¬ 
fer,  but  also  the  temperature  dependence  of  thermoelectric 
properties  is  considered  in  this  paper.  It  is  not  the  full  tempera¬ 
ture  dependence  which  could  only  be  determined  numerically, 
but  it  comes  indirectly  into  play  due  to  averages.  These  averages 
are  constants  and  the  equations  are  solved  by  adapting  constant 
material  properties.  For  given  heat  source  and  heat  sink  tem¬ 
peratures,  Th  and  TL,  the  junction  temperatures  of  thermoelec¬ 
tric  element,  T \  and  Tc,  are  unknown,  thus  the  thermoelectric 
properties,  ah,  ac,k,  a  and  JZ ,  are  unknown.  Alteration  method 
is  adopted  herein  to  determine  the  junction  temperatures.  For 
given  initial  values  of  Th  and  Tc  (Th  =  TH,  Tc  =  TL,  for  exam¬ 
ple),  ah,  ac,  k,  a  and  JZ  can  be  calculated  by  the  fitting  equa¬ 
tion  of  the  thermoelectric  material  (Eqs.  (41)-(43)  hereinafter). 
Then  Th  and  Tc  can  be  calculated.  The  alteration  process  will  be 
repeated  until  required  precision  is  obtained. 

3.  Numerical  example 

A  practical  thermoelectric  generator  in  engineering  is 
adopted  herein  to  analyze  the  effects  of  heat  transfer  laws  on 
the  maximum  power  output  and  maximum  efficiency.  In  the 
numerical  analysis  and  optimization,  the  number  and  size  of  the 
thermoelectric  elements  are  set  as  N  =  127,  A  =  1  x  1  mm2, 
and  L  =  2  mm  [65].  The  heat  reservoir  temperatures  are  set 
as  Th  =  450  K  and  TL  =  300  K.  The  physical  properties  of  the 
commercially  available  material  (Bi2Te3)  by  Melcor  are  shown 
as  follows  [65]: 


L”  =  1  -  (1  -  /)/“ 1  (Tcn  -  T”)m  (T"  -  Th")-m 

+  Xp(Kpf  (T/j  -  T^)m  -  N(ah  ITh  +  I((Th  -  Tc) 
-0.5I2R-0.5JiI(Th  -Tc))) 

+  Mkt(  1  -/)(TC"  -  TnL)m  ~  N(acITc 
+  K(Th-Tc)+0.5I2R  +  0.5JU(Th-Tc))),  (33) 

where  and  0^  are  the  Lagrangian  multipliers,  which  are 
constants  to  be  determined.  The  optimal  ratio  of  thermal 
conductance  allocation,/^,  optimal  electrical  current,  /^,  and  the 
corresponding  hot  junction  temperature,  Th,  and  cold  junction 
temperature,  Tc,  should  satisfy  the  Euler-Lagrange  equations: 


3  L«  dU1  dU1  dU1 

-  =0,  -  =0,  -  =0,  -  =0. 

dTh  3  Tc  3/  3/ 

Eq.  (34)  can  be  calculated  as: 

(1  -f)f-'DcmnT”-'Di'  (T”  -  T £)_1 

-X.i-KjfmnT^Dh  (T”  -  T,")-1 

-  N(ah  I  +  K  -  0.5 JZI))  -  </>nN(K  +  0.5/1/)  =  0, 

(1  -  f)f-'mnT?-'DcD; C  ( Tnc  -  T/1)-1 
+  X^Ni-K  +  0.5JZI)  ~  4>p(N(cic  I  -K  -  0.5JZI) 


(34) 


(35) 


ap  =  (22  224.0  +  930.6T  -  0.9905T2)10~9  VK“\  (41) 

p  =  (5112.0  +  163.4T  +  0.6279T2)10-1°  Om,  (42) 

k  =  (62  605.0-  277.71  +  0.4131T2)10“4  Wm_1  K-1,  (43) 

where  ap,  p  and  k  are  the  Seebeck  coefficient,  electrical 
resistivity  and  thermal  conductivity.  The  Thomson  coefficient 
is  given  by  the  second  Kelvin  relationship  [1-5] 

da  , 

/jl  =  T  —  .  44 

dT 

Figures  2-4  show  the  effect  of  heat  transfer  law  on 
temperature  difference  (AT  =  Th  —  Tc)  between  the  hot  and 
cold  junctions,  power  output  P  and  efficiency  r]  versus  working 
electrical  current  /,  respectively.  Figure  5  shows  the  effect  of 
heat  transfer  law  on  power  output  P  versus  efficiency  rj.  In 
the  calculations,  the  ratio  of  thermal  conductance  allocation 
is  set  as  /  =  0.5.  It  can  be  seen  that  the  heat  transfer  law 
does  affect  the  performance  of  the  irreversible  thermoelectric 
generator.  According  to  the  order  of  temperature  difference, 
they  are  ranked  as  follows:  Special  complex  heat  transfer  law, 
Dulong-Petit  heat  transfer  law,  radiative  heat  transfer  law, 
linear  phenomenological  heat  transfer  law,  and  Newtonian  heat 
transfer.  The  effects  of  heat  transfer  law  on  the  characteristic  of 
power  output  and  efficiency  versus  working  electrical  current 
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Table  1 :  Optimization  results  of  the  optimal  variables  and  the  corresponding  optimal  performance  with  different  heat  transfer  law. 


Transfer  law 

fp 

Ip 

Pmax  (W) 

rip(  10-2) 

fr, 

u 

??max 

MW) 

Newtonian  (n  =  1 ,  m  =  1) 

0.5268 

0.6014 

2.9122 

0.0656 

0.5287 

0.5111 

0.0670 

2.8465 

Linear  phenomenological  (n  =  —  1,  m  =  1) 

0.6201 

0.5996 

2.8982 

0.0655 

0.6221 

0.5099 

0.0669 

2.8333 

Radiative  (n  =  4,  m  =  1) 

0.3873 

0.5967 

2.8769 

0.0653 

0.3478 

0.5061 

0.0667 

2.8092 

Dulong-Petit  (n  =  1,  m  =  1.25) 

0.5277 

0.5730 

2.6827 

0.0639 

0.5313 

0.4886 

0.0652 

2.6246 

Special  complex  (n  =  4,  m  =  1.25) 

0.3837 

0.5651 

2.6222 

0.0634 

0.3857 

0.4818 

0.0647 

2.5653 

Figure  2 :  Effect  of  heat  transfer  law  on  temperature  difference  versus  working 
electrical  current. 


Figure  3 :  Effect  of  heat  transfer  law  on  power  output  versus  working  electrical 
current. 

curves  are  analogous  with  the  effect  of  heat  transfer  law  on 
the  characteristic  of  temperature  difference  versus  working 
electrical  current.  Among  the  several  heat  transfer  laws,  the 
effect  of  Dulong-Petit  heat  transfer  and  special  complex 
heat  transfer  law  are  different  from  other  heat  transfer  laws 
strikingly. 

The  optimization  results  of  the  optimal  variables  and 
the  corresponding  optimal  performance  with  different  heat 
transfer  laws  are  listed  in  Table  1.  Figure  6  shows  the  optimal 
working  electrical  current  IP  and  optimal  ratio  fP  of  thermal 
conductance  allocation  corresponding  to  the  maximum  power 
output  Pmax  with  different  heat  transfer  laws.  Figure  7  shows  the 
maximum  power  output  Pmax  and  the  corresponding  efficiency 
rjP  with  different  heat  transfer  laws.  From  the  table  and 
figures,  the  changing  features  of  optimal  variables  and  optimal 
performance  can  be  concluded.  According  to  the  order  of  the 
optimal  ratio  fP  of  thermal  conductance  allocation,  they  are 


1(A) 


Figure  4:  Effect  of  heat  transfer  law  on  efficiency  versus  working  electrical 
current. 


Figure  5 :  Effect  of  heat  transfer  law  on  power  output  versus  efficiency. 


ranked  as  follows:  special  complex  heat  transfer  law,  radiative 
heat  transfer  law,  Newtonian  heat  transfer  law,  Dulong-Petit 
heat  transfer,  and  linear  phenomenological  heat  transfer  law. 
According  to  the  order  of  the  optimal  electrical  current  IPt 
they  are  ranked  as  follows:  special  complex  heat  transfer 
law,  Dulong-Petit  heat  transfer,  radiative  heat  transfer  law, 
linear  phenomenological  heat  transfer  law  and  Newtonian  heat 
transfer  law.  According  to  the  order  of  the  maximum  power 
output  Pmax,  they  are  ranked  as  follows:  special  complex  heat 
transfer  law,  Dulong-Petit  heat  transfer,  radiative  heat  transfer 
law,  linear  phenomenological  heat  transfer  law  and  Newtonian 
heat  transfer  law.  The  order  of  rjP  is  same  to  the  order  of  Pmax. 
Figure  8  shows  the  optimal  working  electrical  current 
and  optimal  ratio  of  thermal  conductance  allocation 
corresponding  to  the  maximum  efficiency  ^max  with  different 
heat  transfer  laws.  Figure  9  shows  the  maximum  efficiency  r]max 
and  the  corresponding  power  output  P n  with  different  heat 
transfer  laws.  According  to  the  order  of  the  optimal  ratio  of 
thermal  conductance  allocation/^,  they  are  ranked  as  follows: 
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Figure  6:  Optimal  working  electrical  current  and  optimal  ratio  of  thermal 
conductance  allocation. 
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Figure  7:  Maximum  power  output  and  the  corresponding  efficiency  with 
different  heat  transfer  law. 


radiative  heat  transfer  law,  special  complex  heat  transfer  law, 
Newtonian  heat  transfer  law,  Dulong-Petit  heat  transfer,  and 
linear  phenomenological  heat  transfer  law.  According  to  the 
order  of  the  optimal  electrical  current  J^,  they  are  ranked  as 
follows:  special  complex  heat  transfer  law,  Dulong-Petit  heat 
transfer,  radiative  heat  transfer  law,  linear  phenomenological 
heat  transfer  law  and  Newtonian  heat  transfer  law.  According 
to  the  order  of  the  maximum  efficiency  r]maxi  they  are  ranked  as 
follows:  special  complex  heat  transfer  law,  Dulong-Petit  heat 
transfer,  radiative  heat  transfer  law,  linear  phenomenological 
heat  transfer  law  and  Newtonian  heat  transfer  law.  The  order 
of  Prj  is  the  same  as  the  order  of  r]max. 

Since  each  heat  transfer  process  has  a  corresponding  mecha¬ 
nism,  one  could  not  simply  make  a  comparison  just  considering 
the  effect  of  temperature  (the  power  of  temperature  or  tem¬ 
perature  difference).  For  example,  the  I<H  and  I<L  in  the  paper 
represent  heat  transfer  coefficient.  To  convective  heat  transfer 
process,  I<H  and  I<L  represent  convective  heat  transfer  coeffi¬ 
cient;  to  radiation  process,  I<H  and  I<L  represent  a  coefficient 
(involves  Stefan-Boltzmann  constant)  whose  value  is  much 
smaller  than  the  former.  In  the  analysis  mentioned  above,  the 
effect  of  heat  transfer  coefficient  has  not  been  taken  into  ac¬ 
count.  Some  results  only  hold  for  the  numerical  values  adopted 
herein. 


Figure  8:  Optimal  working  electrical  current  and  optimal  ratio  of  thermal 
conductance  allocation  corresponding  to  the  maximum  efficiency  with 
different  heat  transfer  law. 


Figure  9:  Maximum  efficiency  and  the  corresponding  power  output  with 
different  heat  transfer  law. 


of  irreversible  thermoelectric  generator  with  generalized  heat 
transfer  law  Q  oc  A(Tn)m  is  established.  A  class  of  heat  transfer 
laws  is  included  and  multi-irreversibilities  are  considered  in 
the  model.  Applying  the  model  to  a  practical  example  in 
engineering,  it  is  found  that  the  external  heat  transfer  law 
does  affect  the  characteristics  of  the  thermoelectric  device.  It 
is  proved  that  there  are  optimal  working  electrical  currents  and 
optimal  ratio  of  thermal  conductance  allocations  corresponding 
to  the  maximum  power  output  and  maximum  efficiency  when 
the  heat  transfer  law  is  nonlinear.  The  changing  features  of  the 
maximum  power,  maximum  efficiency  and  the  corresponding 
design  variables  are  analyzed  in  detail. 

As  the  heat  exchangers  for  thermoelectric  device  are  various, 
the  heat  transfer  laws  are  various  and  different  from  each 
other.  One  can  determine  the  exponents  n  and  m  in  Q  oc 
A(Tn)m  by  experiment  or  from  empirical  formula.  Combining 
the  heat  transfer  law  and  the  results  obtained  herein  can 
offer  principles  for  the  power  and  efficiency  optimizations  of 
practical  thermoelectric  generators  at  various  external  heat 
transfer  conditions.  The  model  and  optimization  method  may 
be  applied  to  the  analysis  and  design  of  practical  thermoelectric 
generators. 


4.  Conclusions 
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